function [elements,processInfo] = InitializeSolutionStep...
    (elements,model,processInfo)
ne = model.countElements();
for i=1:ne
    e = model.getElement(i-1);
    f = e.getFace();
    ue = f.getUHatLocal();
    xi = e.getGaussPoint(0);
    B = zeros(3,6);
    B = e.computeB(f,B,xi);
    C = elements(i).constitutiveMatrix;
    if elements(i).isCrack == 1
        %compute epsilon
        H = elements(i).H;
        e = elements(i).e;
        epsilon = B*(ue - H*e);
        %compute sigma
        sigma = C*epsilon;
    else
        sigma = C*B*ue;
    end
    elements(i).stress = sigma;
end
end
